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Abstract 

The Casimir energy is evaluated for massless scalar fields under Dirichlet or Neumann boundary 
conditions, and for the electromagnetic field with perfect conductor boundary conditions on one and 
two infinite parallel plates moving by uniform proper acceleration through the Fulling Rindler vacuum 
in an arbitrary number of spacetime dimension. For the geometry of a single plate both regions of 
the right Rindler wedge, (i) on the right (RR region) and (ii) on the left (RL region) of the plate are 
considered. The zeta function technique is used, in combination with contour integral representations. 
The Casimir energies for separate RR and RL regions contain pole and finite contributions. For an 
infinitely thin plate taking RR and RL regions together, in odd spatial dimensions the pole parts cancel 
and the Casimir energy for the whole Rindler wedge is finite. In d = 3 spatial dimensions the total 
Casimir energy for a single plate is negative for Dirichlet scalar and positive for Neumann scalar and 
the electromagnetic field. The total Casimir energy for two plates geometry is presented in the form 
of a sum of the Casimir energies for separate plates plus an additional interference term. The latter 
is negative for all values of the plates separation for both Dirichlet and Neumann scalars, and for the 
electromagnetic field. 

PACS number(s): 03.70.+k, ll.10.Kk 



1 Introduction 



The Casimir effect is a phenomenon common to all systems characterized by fluctuating quantities on 
which external boundary conditions are imposed. It may have important implications on all scales, from 
cosmological to subnuclear. The imposition of boundary conditions on a quantum field leads to the modi- 
fication of the spectrum for the zero-point fluctuations and results in the shift in the vacuum expectation 
values for physical quantities such as the energy density and stresses. In particular, the confinement of 
quantum fluctuations causes forces that act on constraining boundaries. The particular features of the 
resulting vacuum forces depend on the nature of the quantum field, the type of spacetime manifold, the 
boundary geometries and the specific boundary conditions imposed on the field. Since the original work 
by Casimir in 1948 [1] many theoretical and experimental works have been done on this problem (see, 
e.g., [2, 3, 4, 5, 6, 7, 8] and references therein). Many different approaches have been used: mode summa- 
tion method with combination of the zeta function regularization technique, Green function formalism, 
multiple scattering expansions, heat-kernel series, etc. An interesting topic in the investigations of the 
Casimir effect is the dependence of the vacuum characteristics on the type of the vacuum. It is well 
known that the uniqueness of the vacuum state is lost when we work within the framework of quantum 
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field theory in a general curved spacetime or in non-inertial frames. In particular, the use of general co- 
ordinate transformations in quantum field theory in flat spacetime leads to an infinite number of unitary 
inequivalent representations of the commutation relations. Different inequivalent representations will in 
general give rise to different pictures with different physical implications, in particular to different vacuum 
states. For instance, the vacuum state for an uniformly accelerated observer, the Fulling-Rindler vacuum 
[9, 10, 11, 12], turns out to be inequivalent to that for an inertial observer, the familiar Minkowski vacuum 
(for a mathematical discussion by means of a normal mode analysis sec Rcf. [13]). Quantum field theory 
in accelerated systems contains many special features produced by a gravitational field. This fact allows 
one to avoid some of the difficulties entailed by renormalization in a curved spacetime. In particular, the 
near horizon geometry of most black holes is well approximated by Rindler and a better understanding 
of physical effects in this background could serve as a handle to deal with more complicated geometries 
like Schwarzschild. The Rindler geometry shares most of the qualitative features of black holes and is 
simple enough to allow detailed analysis. Another motivation for the investigation of quantum effects in 
the Rindler space is related to the fact that this space is conformally related to the de Sitter space and 
to the Robertson- Walker space with negative spatial curvature. As a result the expectation values of the 
energy-momentum tensor for a conformally invariant field and for corresponding conformally transformed 
boundaries on the de Sitter and Robertson- Walker backgrounds can be derived from the corresponding 
Rindler counterpart by the standard transformation (see, for instance, [14]). 

The problem of vacuum polarization brought about by the presence of an infinite plane boundary 
moving with uniform acceleration through the Fulling-Rindler vacuum was investigated by Candelas and 
Deutsch [15] for the conformally coupled 4D Dirichlet and Neumann massless scalar and electromagnetic 
fields. In this paper only the region of the right Rindler wedge to the right of the barrier is considered. In 
Ref. [16] we have investigated the Wightman function and the vacuum expectation values of the energy 
momentum-tensor for the massive scalar field with general curvature coupling parameter, satisfying the 
Robin boundary conditions on the infinite plane in an arbitrary number of spacetime dimensions and for 
the electromagnetic field. Unlike Ref. [15] we have considered both regions, including the one between 
the barrier and Rindler horizon. The vacuum expectation values of the energy-momentum tensors for 
scalar and electromagnetic fields for the geometry of two parallel plates moving by uniform acceleration 
are investigated in Ref. [17]. In particular, the vacuum forces acting on the boundaries are evaluated. 
They are presented as a sum of the interaction and self-action parts. The interaction forces between 
the plates are always attractive for both scalar and electromagnetic cases. The self-action forces contain 
well-known surface divergencies and needs further regularization. In Refs. [16, 17] the mode summation 
method is used in combination with the generalized Abel-Plana summation formula [18]. This allowed us 
to present the vacuum expectation values in the terms of the purely Rindler and boundary parts. Due to 
the well known non-integrable surface divergences in the boundary parts, the total Casimir energy cannot 
be obtained by direct integration of the vacuum energy density and needs an additional regularization. 
Many regularization techniques are available nowadays and, depending on the specific physical problem 
under consideration, one of them may be more suitable than the others. In particular, the generalized zeta 
function method [19, 20] is in general very powerful to give physical meaning to the divergent quantities. 
There are several examples of the application of this method to the evaluation of the Casimir effect (see, 
for instance, [7, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]). In this paper, by using the zeta function 
technique, we consider the Casimir energy for the geometries of a single and two parallel plates, moving 
by uniform proper acceleration through the Fulling-Rindler vacuum. 

The paper is organized as follows. In Sec. 2 the Casimir energy is evaluated for the Dirichlet scalar 
in the case of a single plate. For the both RR and RL regions we construct integral representations of 
the related zeta functions and analytically continue them to the physical region. The geometry of two 
plates in the case of the Dirichlet scalar is investigated in Sec. 3. We show that the corresponding zeta 
function is a sum of the zeta functions for the separate plates plus an additional interference term which 
is finite in the physical region. Similar problems for the Neumann scalar are investigated in Sec. 4 and 
Sec. 5, respectively. Section 6 considers the Casimir energy for the electromagnetic field assuming that 
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the plates are perfect conductors. Section 7 concludes the main results of the paper. In Appendix A the 
case d = 1 is considered separately. Appendix B discusses the relation between the calculated energies 
and the energies measured by an uniformly accelerated observer. 

2 Casimir energy for a single Dirichlet plate 

Consider a real massless scalar field ip(x) with curvature coupling parameter ( satisfying the field equation 

V M VV + Ci^ = 0, (2.1) 

with R being the scalar curvature for ad + 1-dimensional background spacetime, is the covariant 
derivative operator associated with the corresponding metric tensor g^ v . For minimally and conformally 
coupled scalars one has ( = and ( = (d — l)/4d, respectively. Our main interest in this paper will be 
the Casimir energy in the Rindler spacetime induced by a single and two parallel plates moving with 
uniform proper acceleration when the quantum field is prepared in the Fulling- Rindler vacuum (note the 
difference of our problem from that related to the Unruh effect, where the field is in its Minkowski vacuum 
state). For this problem the background spacetime is flat and in Eq. (2.1) we have R = 0. As a result 
the eigenmodes are independent of the curvature coupling parameter and, hence, the Casimir energy will 
not depend on this parameter. However, the local characteristics of the vacuum such as energy density 
and vacuum stresses depend on the parameter £ [16, 17]. 

In the accelerated frame it is convenient to introduce Rindler coordinates (r, £, x) which are related 
to the Minkowski ones, (t, x 1 ,x) by transformations 

t = £sinhr, x 1 = £coshr, (2-2) 

where x = (x 2 , . . . , x d ) denotes the set of coordinates parallel to the plate. In these coordinates the 
Minkowski line element takes the form 

ds 2 = i 2 dT 2 - di 2 - dx 2 , (2.3) 

and a wordline defined by £,x = const describes an observer with constant proper acceleration 
Rindler time coordinate r is proportional to the proper time along a family of uniformly accelerated 
trajectories which fill the Rindler wedge, with the proportionality constant equal to the acceleration. 
Let {ip a (x), ip* a (x)} be a complete set of positive and negative frequency solutions to the field equation 
(2.1), where a denotes a set of quantum numbers. For the geometry under consideration the metric and 
boundary conditions are static and translational invariant in the hyperplane parallel to the plates. It 
follows from here that the corresponding part of the eigenfunctions can be taken in the standard plane 
wave form: 

ip a = C </)(£) exp [i (kx - lot)] , a = (k,u), k = (fc 2 , ■ ■ ■ , k d ). (2.4) 
The equation for </>(£) is obtained from field equation (2.1) on background of metric (2.3): 

evco + w(o + - k 2 e) m = o, (2.5) 

where the prime denotes a differentiation with respect to the argument, and k = |k|. The linearly 
independent solutions to equation (2.5) are the Bessel modified functions Iiuj(k^) and Ki U! (k^) of the 
imaginary order. The eigenfrequencies are determined from the boundary conditions imposed on the 
field on the bounding surfaces. 
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Figure 1: The (xi,t) plane with the Rindler coordinates. The heavy line £ = £1 represent the trajectory 
of the plate. 

2.1 Vacuum energy in the RR region 

In this section we will consider the vacuum energy for a scalar field satisfying Dirichlet boundary condition 
on a single plate located at £ = £1: 

Vl€=€i = 0. (2.6) 

We will assume that the plate is situated in the right Rindler wedge x l > \t\. The surface £ = £i 
represents the trajectory of the boundary, which therefore has proper accelerations £J~ 1 (see Fig. 1). This 
trajectory divides the right Rindler wedge into two regions with £ > £i and £ < £i- In the following we 
will refer these regions as RR and RL regions, respectively. First let us consider the vacuum energy in 
the RR region. In this region, for a complete set of solutions that are of positive frequency with respect 
to d/dr and bounded as £ — > oo, in Eq. (2.4) one has to take </>(£) = ZQ W (A;£). For the Dirichlet scalar 
the corresponding eigenfrequencies are determined from the equation 

KUHi) = 0. (2.7) 

The positive roots to this equation arranged in the ascending order we will denote by u> = u>iDn(k£i), 
n = 1,2, . . ., oJiDn < wiDn+i- To distinguish in the notations we will use indices D and N to denote 
Dirichlet and Neumann boundary conditions respectively. Similarly the quantities for the RR and RL 
regions will be denoted by indices R and L, respectively. The vacuum energy per unit surface of the plate 
in the region £ > £i is given by formula 

= o / otE^(^i) = -A / dxx d - 2 ^Dn(x), (2.8) 

ZJ V*) n=l £l J ° n=l 

where we have integrated over angular coordinates, and 

Bd = - ,J - - ■ (2.9) 
(4vr) — r(^fi) 

Here and below we will assume that d > 1. The case d = 1 is considered in Appendix A. As it stands, 
the right-hand side of equation (2.8) clearly diverges and needs some regularization. We regularize it by 
defining the function 

/•oo 

Zdk(s) = / dxx d - 2 C£l(8,x) (2.10) 
Jo 
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where we have introduced the partial zeta function related to the eigenfrequences defined by Eq. (2.7): 



71=1 



The Casimir energy is expressed as 



-^Z DK {s) 
si 



-i- 



(2.11) 



(2.12) 



In accordance with this formula, the computation of the Casimir energy requires the analytic continuation 
of the zeta function to the value s = — 1. 

The starting point of our consideration is the representation of the partial zeta function for the 
corresponding modes in term of contour integral (for a similar treatment of the zeta function as a contour 
integral see Refs. [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]): 



(2.13) 



where C is a closed counterclockwise contour in the complex z plane enclosing all zeros 0J\Dni, x )- We 
assume that this contour is made of a large semicircle (with radius tending to infinity) centered at the 
origin and placed to its right, plus a straight part overlapping the imaginary axis. When the radius of 
the semicircle tends to infinity the corresponding contribution into Q D '(s,x) vanishes for Res > 1. After 
parametrizing the integrals over imaginary axis we arrive at the expression 

AR), \ 1 . f°° 7 q d . . 

C 1D ( s > x ) = ~ sm y J dzz g^ lnK z{x)- 



(2.14) 



For small z, z — > 0, one has K z (x) m K (x) + (d 2 K z (x)/dz 2 ) z= o(z 2 /2), and the integral in (2.14) 
converges at the lower limit for Res < 2. Hence, this integral representation of the zeta function is valid 
for 1 < Res < 2. To evaluate the Casimir energy we need the zeta function at s = — 1 and so we have to 
do an analytic continuation of Eq. (2.14). To this aim we employ the uniform asymptotic expansion of 
the modified Bessel function K z (x) for large values of the order [33]: 



K z {x) 



n e -ZT)(x/z) 

2(x 2 + z 2 ) 1/ * 



Ki D \x), tff>(z) ~jT(-l) 



ui(t) 



1=0 



where 



rj(x) = y/l + x 2 + In ■ 



ui{t) 



(x 2 + z 2 ) 112 ' 
ui(t) 



(2.15) 



(2.16) 



Vx 2 + z 2 l + VTT^ 2 ' 

and the expressions for the functions ui(t) are given in [33]. From these expressions it follows that the 
coefficients ui(t) have the structure 



2m 



(2.17) 



m=0 



with the numerical coefficients ui m . From the recurrence relations for the polynomials ui(t) (see [33]), 
the following recurrence formulae can be obtained for the coefficients u\ m : 



1 



Ul+l,m 



2m + l + 



4(2m + / + 1) 



where m = 0, 1, ...,/ + 1, and = u\j + i 
zeta function can be presented in the form 



1 

0, 1i 



2m + / - 2 + 



(2.18) 



4(2m + / + 1)_ 
1. Now making use the expansion (2.15), the 



cS(«,x)=c 1 ( s o) ( S ,x)+c;? i) ( 5 ,x) 



(2.19) 
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where 



rz?m 1 its f°° d fir er z ^ x l z ) 

<$"(», x) = ±sin^°°,fez-«^lnirP(x). (2.21) 
Under the conditions 1 < Res < 2, the term (2.20) can be easily evaluated to give 



Jm), \ (x/2) 1 s f s 1 \ . its x- 

Qr, (s,x) = — ; -B\l — s, sin , (2.22) 

SlD v ' ' vr(l-s) V 2 / 2 4 



with the beta function B(x,y). Now using the standard dimensional regularization result that the 
renormalized value of the integrals of the type Jq°° dxx@ is equal to zero (see, e.g., [34]), we conclude 
that the contribution of the term Qp(s,x) into Eq. (2.10) vanishes. This can be seen by another 
way, considering the case of a scalar field with nonzero mass m and taking the limit m — ► after the 
evaluation of the corresponding integrals (for this trick in the calculations of the Casimir energy see, for 
instance, Refs. [35, 36]). For the massive case in Eq. (2.4) one has </>(£) = K^^s/k 2 + m 2 ) and the 
corresponding formulae for Q D '(s,x) are obtained from those given above in this section by replacement 

x — > x 2 + ?rt 2 £ 2 . With this replacement the integral corresponding to the contribution of Cid°^ ^ n *° 
Eq. (2.10) can be easily evaluated in terms of gamma function and vanishes in the limit m — ► for 
Re s < d — 1. For this reason in the following we will concentrate on the contribution of the second term 
on the right of Eq. (2.19). The corresponding expression for Zdk{s) from Eq. (2.10) takes the form 

1 POO POO Q 

Z DK (s) = -sin^f / dxx d - 2 dz z- 3 -^- In K<i D \x). (2.23) 
tt 2 J J dz 

To deal with one infinite limit integral in the following it will be convenient to introduce polar coordinates 
on the plane (z, x): 

z = rcos9, x = rs'm9, (2.24) 



In these coordinates one has the following expressions 

t = co S 9, z V (x/z)=rg(9), K™ (x) ~ |^(-i) ^( c ° s ^ . ( 2 .25) 

1=0 

where we have introduced the notation 



sin B 

g(0) = l + cos 9 In . (2.26) 

yw l + cos(9 y ' 

The integral representation (2.23) is well suited for the analytic continuation in s. Following the usual 
procedure applied in the analogous calculations (see, for instance, [7] and references therein), we subtract 
and add to the integrand in Eq. (2.23) M leading terms of the corresponding asymptotic expansion and 
exactly integrate the asymptotic part. For this let us define the coefficients Ui(cos9) in accordance with 

E° , w Ui( cos 9) 
{ - l) Trrw^- (2 - 27) 

l=i (l + r 2 ) 7 

Here on the right we have expanded over (1 + r z ) to avoid convergence problems in the integral over 
r at the lower limit r = 0. Hence, Eq. (2.23) may be split into the following pieces 

ZDK{8) = Z<${a) + zg> K {8), (2.28) 
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with 

M 



r,(as), \ 1 7TS f°° , d 2 f°° , J A, 1W Ul(cOs6) 

ZU(S) = -sm T y o cfa*--^ *,-_£(-l)'-iL^, (2.29) 



7 (i) 



(.) = Is^^^^^^^^^w-g^,.^). (,30, 



The number of terms M in these formulae is determined by the spatial dimension of the problem under 
consideration. For M > d the expression Z^ K (s) is finite at s = — l and, hence, for our aim it is 
sufficient to subtract M = d asymptotic terms. For larger M the convergence is stronger. This means 
that with M > dwe can directly put s = — l in the expression for Z^ K (s) in Eq. (2.30) and perform the 
integral numerically. For numerical evaluation it is useful to integrate by parts the z-integral and then 
to introduce polar coordinates in the (z, x)-plane. This yields 

M 



= 1 Hdrr^ r /2 d9 sm^e In ( ^% rasJ (r sinfl) ) -£(-1)1^5*) 
n Jo Jo \ V vr J fit {l + r 2 ) l,A 

(2.31) 

Now we turn to the asymptotic part of Zdk(s) given by expression (2.29). For this term the analytic 
continuation can be done analytically because its structure is simple. To evaluate this part we integrate 
over z by parts and note that the functions Ui(t) are polynomials in t: 

i 

Um = £ U lm t 2m , (2.32) 

m=0 

where the numerical coefficients Ui m are related to the coefficients in Eq. (2.17) by expansion (2.27). 
The first five functions are 

UAt) = --—t\ u 2 {t) = — --t 2 + —t\ (2.33a) 

w 8 24 ' w 16 8 16 

„ , , 49 1793 2 221 4 1105 fi 

U 3 (t) = t 2 H t 4 1 6 , (2.33b) 

SK ' 384 1920 128 1152 ' V ; 

„ , , 21 83 2 551 4 339 fi 565 8 

UAt) = 1 2 + t 4 f + t 8 , 2.33c 

v ; 128 32 64 32 128 ' V ' 

1813 297649 2 198755 4 136907 fi 82825 s 82825 10 

U 5 (t) = 1 2 H 1 4 1 6 H t 8 1 10 . (2.33d) 

K ' 5120 35840 4608 1536 1024 3072 V ; 

Taking into account Eq. (2.32) and introducing polar coordinates, for the asymptotic part one finds 

Z<#(«) = £^£(-1)' E Ui m / dr / / d6 sin"" 2 cos 2 — 1 9. (2.34) 

* 1 l=1 m=0 Jo (l + r 2 )' Jo 

Evaluating the integrals by using the standard formulae (see, for instance, [37]), we arrive at the expression 
<»M = ^i-fB-D' E U lm B (IZIZI, ! + ; - j " + ' ) L - • , (, 35 ) 

(=1 m=0 ^ / V / 

where the pole contributions are given explicitly in terms of beta function. In the sum over I, the terms 
with even d — I > have simple poles at s = — 1 coming from the first beta function. Introducing a new 
summation variable p = (d — l)/2, the corresponding residue can be easily found by using the standard 
formula for the gamma function: 

( l\d d ~ 2 P R ( m i 1 d-l\ 

p=0 m=0 v 2 / 
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Table 1: Pole and finite parts of the Casimir energy for a single Dirichlet plate. 



d 


cd-l E (H) 
SI ^IDp 


SI ^lDf 


SI ^lDp 


Si ^lDf 


l£; iD/ 


2 
3 

4 


1 


0.00185 
-0.000620 
0.000225 


1 


-0.00269 
-0.000695 
-0.000235 


-0.000836 
-0.00131 
-0.00000996 


512tt(s+1) 
1 


512tt(s+1) 


1260tt 2 (s+1) 


1260tt 2 (s+1) 
37 


262144vr 2 (s+l) 


262144tt 2 (s+1) 



Pd 



where 

(d-l)/2, 
d/2 - 1 

Hence, Laurent-expanding near s = — 1 we can write 



for odd d 
for even d 



with 



7 (as) _ 



(-1) 
2vrd 



p=0 



d-2p 

-if £ ^ 

m=0 



■ J DK,-1 
S + l 



+ zjs$ J0 + o{8 + i), 



(2.37) 



(2.38) 



2p,m , a , -, ■> 




m + \ 



- 2 



(-1)' y, u lm B 



d I 
2'~ 



B \m + 



1 d-1 



(2.39) 



where i/)(x) = d\aT{x)/dx is the diagamma function and the second sum in the braces of the third line 
is present only for M > d + 1. The first term on the right of Eq. (2.39) with diagamma functions 
comes from the finite part of the Laurent expansion of the summands with even d — I > in Eq. (2.35). 
Gathering all contributions together, near s = — 1 we find 



Zdk(s) 



7 (as) 
J DK,- 



8 + 1 



i + 4^ + Z^(-l)+O( S + l), 



(2.40) 



where the separate terms are defined by formulae (2.31), (2.36), (2.39). Using this result, for the vacuum 
energy induced by a single plate at £ = £i in the region £ > £i one receives 



where for the pole and finite contributions one has 



^IDp + ^lDf 



E 



(R) 
lDp 



B d Z 



(as) 

DK,-1 



E 



(R) 
IDf 



Ba 



7 {as) 
'DKfl 



+4'k-i) 



(2.41) 



(2.42) 



ttHs + iy 1Dt it 1 

The numerical results corresponding to the vacuum energy (2.41) with separate parts (2.42) are presented 
in the first two column of Table 1 for spatial dimensions d = 2, 3, 4. 

In Appendix B we have discussed the relation between vacuum energy (2.41) and the energy measured 
by an uniformly accelerated observer with the proper acceleration g. These energies are connected by 
formula (B.4), where ji is an arbitrary constant with the dimension of mass. Using the numerical results 
given in Table 1, the Casimir energy in the RR region for the massless Dirichlet scalar in d = 3 measured 
by an uniformly accelerated observer is presented in the form 



-0.000620 + 



1 



1260vr 2 



1 



s+l 



+ In 



9 J J 



(2.43) 
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where the logarithmic term is a consequence of the divergence and has to be viewed as a remainder of 
the renormalization process. The discussion for the role of the normalization scale \i in the calculations 
of the Casimir energy can be found in Ref. [21]. The remained pole term in the Casimir energy is a 
characteristic feature for the zeta function regularization method and has been found for many cases of 
boundary geometries. Note that a very powerful tool for studying divergence structure of the vacuum 
energy is the heat kernel expansion (see, for instance, [7, 38] and references therein). In particular, the 
coefficient of the logarithmic term is determined by the corresponding boundary coefficient in the heat 
kernel asymptotic expansion. 



2.2 Vacuum energy in the RL region 

Now we turn to the Casimir energy in the region between a single plate and the Rindler horizon corre- 
sponding to £ = 0. As in previous subsection we will assume that the plate is located at £ = £i and the 
field satisfies Dirichlet boundary condition (2.6) on it. In the next section we show that the corresponding 
vacuum energy can be determined by the formula 

£$(&) = ^k/(')|.=-i, (2.44) 
si 

where we have defined the function 

f 00 

Z DI (s)= / dxx^C^&x), (2.45) 
Jo 

with the partial zeta function 

CS (a, x) = ± sin ^ dz z- s ^- In I z (x). (2.46) 

Here p is a small number which we will put zero later. Unlike to the case (2.14), here we can not directly 
put p = 0, as for values s with convergence at the upper limit, the z-integral will diverge at the lower 
limit. The analytical continuation for (2.46) as a function on s can be done by the way similar to that 
given above for the region on the right of a single plate. First, we note that for the Bessel modified 
function one has the uniform asymptotic expansion [33] 

I,{x) = -}= e ^' 4 /f)(x), /P^-E Mt) , /2 , (2.47) 

where the functions rj(x), t, ui(t) are defined in Eq. (2.16). This allows us to present the partial zeta 
function in the form 

(^(s,x)=({ L D \ S ,x) + ({ L D 1 \s,x), (2.48) 

with 

rrn\ 1 7TS f'°° d pzn(x/z) 

dn\s,x) = -sin— / dzz- s ^-\n — , (2.49) 

t[ L D l \s,x) = Isin^^dzz-^ln/f)^). (2.50) 

For 1 < Res < 2 the z-integral in the expression for (^^\s,x) converges in the limit p — > and in this 
limit the evaluation of the integral gives 

C 1D x) = -^^ B (1 ~ ^ —) «n y - — • (2.51) 
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As in the case of the RR region the contribution of this term into the vacuum energy vanishes. 
Now we consider the part due to Eq. (2.50): 

Z DI (s) = -sm^- f°°dxx d ' 2 f°° dzz~ s ^- In I ( z D) (x). 

7T 2 Jo J p OZ 

By using asymptotic expansion (2.47), we can present it in the form 

Z DI (s) = Z^\s) + Z^(s), 



where 



Z^(s) = -sm- 



o 



M 



1 • ™ r dxx^ rw« a ^ ui{cose) 



1 • ™ f°° , d-2 

— sm — / dxx 
2 J 



d 



M 



Ui(cos9) 



2^/2 



(2.52) 
(2.53) 

(2.54) 
(2.55) 



In the expression for Z^\s) we can directly put p = 0, and after introducing polar coordinates and 
integrating, one finds 



1=1 m=0 



S . TVS 

4tt 2 



M I 



d- s-1 l+s-d+1 



2 ' 2 

At s = — 1 this function has a simple pole with the residue 

y{ as ) _ l -\\d 7 {as) 

,(as) 



s d- 1 



2' 2 



(2.56) 



(2.57) 



where Z D ^_ X is determined by formula (2.36). Now we have the following Laurent expansion 

„(as) 

M*)( a \ _ DI-l , 7 {as) 



where 



7 (as) 
J DI,0 



w^" 1 ) 2.^™ B( d_ p>p + 1) 



(2.58) 



p=0 



m=0 



V'(p+l)-V' 7T +V' "l+o -V ; "T-+o 



+ 



47T 



d-1 M 

E + E)E^ 

.Z=l,d-Z=odd l=d+lj m=0 



d l-d 



2' 2 



1 d-1 
B [m+ -, 

' 2 2 



(2.59) 



As regards to Z^j(s), for M > d the corresponding expression (2.54) is finite at s = — 1 and we can 

find 

C//(cos 0) 



directly put p = 0. After integrating by parts and introducing polar coordinates we find 

,tt/2 r M 



Zg}(-1) = - Hdrr^ 1 d8sm d - 2 9 In (v&r rfl Wj rc0B( j(r sin 0)) - V — 
vr Jo Jo v ' f~f (1 



t (1 + r 2 ) 



2\*/2 



(2.60) 



Collecting all terms for the expansion near s = — 1 one finds 



(2.61) 
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with the separate terms defined by (2.57), (2.59), (2.60). Now the vacuum energy in the region < £ < £i 
is the sum of the pole and finite terms 



^S(d) 



E {L) +E (L) 

^lDp + -^LD/' 



with 



^lDf 



B d 

d-i 



z 



(as) 
DI,0 



+ 4 1 J(-i) 



(2.62) 



(2.63) 



The results of the numerical evaluations for these quantities are given in Table 1 for spatial dimensions 
d = 2, 3, 4. The corresponding vacuum energy measured by an uniformly accelerated observer is 

related to the energy (2.62) by formula (B.4). In particular, using the data from Table 1, for the case of 
spatial dimension d = 3 one finds 



ElDgKW ~ 72 
si 



-0.000695 



1260vr 2 



s + l \g 



(2.64) 



where g is the proper acceleration of the observer. 



2.3 Total Casimir energy for a single Dirichlet plate 

The total Casimir energy for a single Dirichlet plate can be obtained by summing the RR and RL parts 
considered above: 



As separate summands, it contains pole and finite parts, 



Eid{£i) — EiDp + EiDf, 



where 



El Dp 



tit 1 (s + l) 



E 



lDf 



Cl 



d-1 



J DK,0 



r(as) 
"DL0 



,(1) 



,(1), 



(2.65) 
(2.66) 

(2.67) 



Note that in odd spatial dimensions the pole part vanishes due to the cancellation of corresponding RR 
and RL parts and the total Casimir energy is finite. In this case this energy can be presented in the form 



E 



ID 



"SI I 1=1 v / m=0 V 



1 d-1 



2' 2 



d-1 1 



+ I drr 
'o 



de sin d - 2 e 



Mi 



ln(2rI rcose (rsm9)K rcose (rsm9)) - 2^ 



U 2 i(cose) 



l= t (1 + r 2 )' 



(2.68) 



For the convergence of the integral in this formula it is sufficient to take Mi > d/2 — 1. The numerical 
results for the finite part of the total vacuum energy are given in Table 1. For the total Casimir energy 
measured by an uniformly accelerated observer with the proper acceleration g, in odd dimensions the 
divergencies cancel and so do the logarithmic terms. In particular, for d = 3 Dirichlet scalar the total 
Casimir energy 

E 1DS = -^ (2.69) 

SI 

is negative. This means that the resulting vacuum forces tend to accelerate the Dirichlet plate. 
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3 Casimir energy for two parallel Dirichlet plates 



In this section we consider the scalar vacuum in the region between two plates located at £ = £1, £ = £2, 
which therefore have proper accelerations ^j -1 and . The problem geometry is depicted in Fig. 2. The 
scalar field satisfies the Dirichlet boundary conditions on the plates: 

<P\t=ti = Vlf=6 = 0- 

The function 0(£) in Eq. (2.4) satisfying the boundary condition on the plate £ = £2 is in the form 




Figure 2: The problem geometry in the case of two plates. The heavy lines £ = £1 and £ = £2 represent 
the trajectories of the plates. 



<K0 = D iuJ (k£, H 2 ) = IUH 2 )K ILU (H) - hu{H)KUH2). (3.2) 

From the boundary conditions on the plate £ = £1 we find that the corresponding eigenfrequencies are 
roots to the equation 

D itJ (k^,k^) = 0. (3.3) 

This equation has an infinite set of solutions. We will denote the positive roots by u = wc„(fc(i, ^£2); 
and will assume that they are arranged in the ascending order toon < WDn+b For the vacuum energy in 
the region £1 < £ < £2 per unit surface of the plates on has 

En = - 2 J ^pE^(^,^). (3.4) 

This energy can be written as 

/>oo 

E D = B d dkk d -\ D (s,kCi,kC 2 )\s=-i, (3.5) 
J 

with the partial zeta function 

00 

CDM£i,fck) = X>^(fc6,fck), Res>l. (3.6) 

n=l 
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In order to obtain the Casimir energy, one has to find the analytic continuation of Eq. (3.6) to s = — 1. 
This is done by an analytic continuation of an adequate contour integration in the complex plane. An 
immediate consequence of the Cauchy's formula for the residues of a complex function is the expression 



Cd{s. 



(3.7) 



where the contour C is the same as in Eq. (2.13) with the additional semicircle C p of small radius p 
avoiding the origin from the right (in Eq. (3.7) we can directly put p = 0, however this semicircle is 
needed for the next step, see Eq. (3.8) below). Let us denote by C + and C~ the upper and lower halves 
of the contour C except the parts coming from the semicircle C p . The integral on the right of Eq. (3.7) 
can be presented in the form 



CD(s,fc£i,fc6) 



+ 



1 



E 



a=+,— 



dz z~ 



c° 



• 9 , r 
' In T 

dz oaz 



+ 



+ 



2ni 



2m 



E 



dz z~ 



,d_ 

dz 



In 



1 - 



Mi)K iz {K 2 ) 
M2)K lz {K l )\ 



dz z 



Co 



dz 



In 



(3.8) 



where we have used (2.13) to introduce the function (^(s,x). The last integral over C p vanishes in 
the limit p — > for Res < 2 and we will omit it in the following consideration. After parametrizing the 
integrals over the imaginary axis, we arrive at the formula 



CD(s,k£i,k€ 2 ) = 



Ci{?(a,fc£i) + - sin ? 

7T I 



dz z 



1 ITS 

H — sin — 

7T 2 



OO Q 

dzz s — In 

oz 



-°-?-lnI z (k&) 
oz 

h(K2)K z {k^) 



(3.9) 



The last integral on the right of this formula is finite at s = — 1 and vanishes in the limits £i — > 
or £2 — ► 00 • It follows from here that the second term on the right corresponds to the zeta function 
Cid ( s 5 ^£2) for the region on the left of a single plate located at £ = £2- The analytic continuation of this 
function we have done in the previous section. 

As a result the vacuum energy in the region £1 < £ < £2 can be written in the form of the sum 



E D = ^(Ci) + 42(6) + A^(Ci,6), 



where the interference term is given by the formula 



TD 1' OC /*CO J 

AE D {£ 1 ,S 2 ) = — dkk d - 2 dzln 1-- 
n Jo Jo I L . 



h{Ki)K z {ki 2 ) 



M2)K z (kCi)\ " 



(3.10) 



(3.11) 



Note that in the corresponding expression we have partially integrated the z-integral. The function 
D z (k^i,k^ 2 ) is positive for £1 < £2 and the energy (3.11) is always negative, AE£>(£i, £2) < 0. In Fig. 3 
we have presented the dependence of the interference part of the Casimir energy on the ratio £1/^2 f° r 
d = 3. In the case of two plates geometry, to obtain the total Casimir energy, , we need to add to 

the energy in the region between the plates, given by Eq. (3.10), the energies coming from the regions 
£ < £1 and £ > £2- As a result one receives 



E 



(tot) 
D 



£id(£i) + Ad (£2) + A£x>(£i,£ 2 ), 



(3.12) 



where the interference part is given by formula (3.11). Note that in the limit £2 — > 00 and for a fixed £1, 
the last two terms on the right vanish and we recover the result for a single plate located at £ = £1. For 
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0.2 0.4 0. 6 0.B 

Figure 3: Interference part of the Casimir energy in the region between two Dirichlet plates, 
£2 _1 A£x>(6-, £2)) as a function on the ratio £i/£2 for d = 3. 

d = 3 Dirichlet scalar all terms on the right of formula (3.12) are negative and, hence, the total Casimir 
energy for two parallel plates is negative. 

In the limit £1 — > £2, expression (3.11) is divergent and for small values of £2/^1 — 1 the main 
contribution comes from large values of z. Replacing the order of integrations and introducing a new 
integration variable x = k/z, we can replace the Bessel modified functions by their uniform asymptotic 
expansions for large values of the order. After evaluating the integrals, to the leading order one receives 

A^(d,6) ~ -(47r)-^0^ + l)r (^) (6 5 6)d , (3.13) 

where Cr( x ) 1S the Riemann zeta function. Note that the energy AEo is related to the corresponding 
quantity AEpg measured by an uniformly accelerated observer by formula AEr, g = gAEp, where g is 
the proper acceleration of the observer (see Appendix B). In the limit £i,£2 —> °°> £2 — £1 = const, 
the parts in Eq. (3.12) corresponding to the contributions from the single plates vanish, -Eii)(£j) ~* 0> 
i = 1,2, and the interference part remains only. In this case, as it follows from Eq. (3.13), the energy 
measured by an observer with the acceleration g — > such that g£i — > 1, coincides with the standard 
Casimir energy for two parallel plates in d+ 1-dimensional Minkowski spacetime [39]. 

As it has been shown in Ref. [17], the vacuum forces acting on the boundaries can be presented as 
a sum of the self-action and interaction forces. It can be easily checked that the interaction forces are 
related to the energy (3.11) by the standard thermodynamical relations 

(7) 

where P^f int y J = 1,2 are the vacuum effective pressures on the plate at £ = £,- induced by the presence 
of the second plate (interaction force per unit surface). 
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4 Casimir energy for a single plate with Neumann boundary condition 
4.1 RR region 

In this section we will consider the vacuum energy for a scalar field satisfying Neumann boundary con- 
dition on a single plate located at £ = £1: 



dip 



e=a = o. (4.1) 



In the region £ > £i, the corresponding eigenfrequencies are roots to the equation 

KUKi) = 0, (4.2) 

where the prime denotes the derivative with respect to the argument of the function. The positive 
roots to this equation arranged in the ascending order we will denote by oj = ooiNn(k£i), n = 1,2, . . ., 
UiNn < wijvn+i- The vacuum energy per unit surface of the plate in the region £ > £i is given by formula 

E m(Zi) = 7/1 / dxx d - 2 J2^Nn(x). (4.3) 

£l J ° n=l 

To regularize this energy we consider the function 

Z NK (s) = / dxx d - 2 cS(s,x), (4.4) 

J 

with the partial zeta function related to the corresponding eigenfrequences: 

oo 

71=1 

The Casimir energy is obtained by the analytic continuation to s = — 1: 

E < Sj{ti) = 4± l Z NK {s)\ s= _ 1 . (4.6) 
si 

This analytic continuation procedure is similar to that presented above for the Dirichlet case. First, we 
present the partial zeta function in terms of a contour integral 

C$(*,z) = ^-J c dzz- s ^-\nK' iz {x), (4.7) 

with the same contour C as in Eq. (2.11). When the radius of the large semicircle of this contour tends 

( Ft} 

to infinity the corresponding contribution into Q N '(s,x) vanishes for Res > 1. After parametrizing the 
integrals over imaginary axis we obtain the following integral representation 

(s, x) = X - sin ^ ^ dz z- s A l n K ' z (x) , (4.8) 

valid for 1 < Res < 2. For the analytic continuation of Eq. (4.8) to s = —1 we use the uniform asymptotic 
expansion of the function K' z (x) for large values of the order. We will write this expansion in the form 

— ( 2 _l_ 2^1/4 oo _ , s 

T 1 ] e-^Kf\x), *W(s) -£(-!)' - Vlit) (4.9) 
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where 



vi(t) 



(4.10) 



and the expressions of the functions vi(t) can be found in Ref. [33]. From these expressions we can see 
that the coefficients vi(t) are polynomials in t: 



vi(t) = ^2 vi m t 



2m 



(4.11) 



m=0 



From the formulae for the functions vi(t) (see [33]), the coefficients vi m can be expressed via the coefficients 
U[ m in Eq. (2.17) by the relations 



1 

--jU lm 



2m + l + l 



1 



4(2m + / + 1) 



1 



2m + 1-1 



4(2m + / + l) 



(4.12) 



where m = 0, 1, . . . , I + 1, and, as has been noted before, ui_\ = li/^+i = 0, u 00 = 1- Using expansion 
(4.9), the zeta function may be split into two pieces 



4W) 



(ra), 



with 



CijV \ s , x ) 



1 


ITS 


f'GG 


- s d_ 




sin — 






7T 


2 


Jo 


dz 


1 


7TS 


/•CO 


-± 




sin — 


/ 




7T 


2 


Jo 


dz 



-ZT](x/z) 



2 x 
ln^f'(x). 



Under the conditions 1 < Res < 2, evaluating the integral in Eq. (4.14) one finds 

\l-s 



C»(^) = <x/2 



„ , , S — 1 \ . 7TS X" 

, is 1 — s, sm 1 

vr(l-s) V 2 / 2 4 



(4.13) 

(4.14) 
(4.15) 

(4.16) 



By the arguments similar to those for the Dirichlet case, we conclude that the contribution of the term 
C[^\s,x) into Eq. (4.4) vanishes. The expression with the second term on the right of Eq. (4.13) takes 
the form 

1 f oo /*oo .CI 

Z NK (s) = -sin— / dxx d ~ 2 dz z~ s hi Kf\x). (4.17) 

7T 2 Jo JO ^ z 

Now we subtract and add to the integrand in Eq. (4.17) the corresponding asymptotic expression and 
exactly integrate the asymptotic part. For this we introduce polar coordinates in the (z, x) plane and 
the coefficients VJ(cos#) in accordance with 



\i=o r J i=i (1 + r 2 ) 7 



This allows us to present (4.17) in the form of the sum 



Z NK (s) = Z^(s) + Z^ K (s), 



(i) 



where 



1 ITS 

— sm — 

7T 2 



7 (l)N 
J NK 



1 . ITS [°° 

— sm — / axx 
* 2 J 



/ dxx d ~ 2 
Jo Jo 

d-2 f 

Jo 



M 

-^J 1 



t d_ 

dz 



i=i 



Vi(cose) 



d 

dzZ ~ S dz 



(l + r 2)« 
M 

ln^f)(*)-J>l) 



*5(cosfl) 
(l + r 2 ) ,/2 



(4.18) 

(4.19) 

(4.20) 
(4.21) 
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For M > d the expression zff K (s) is finite at s = — 1 and we can perform the integral numerically. 
For numerical evaluation it is useful to integrate by parts the z-integral and to perform to the polar 
coordinates. This yields to the final expression 



- / drr d - x / 
71" Jo Jo 



d9 sin 



d-1 , 



(4.22) 



Now let us consider the asymptotic part of the function Zjyx( s ) given by expression (4.20). To 
evaluate this part note that the functions Vi(t) have the structure 

i 

Vi{t) = Vi™t 2m , (4-23) 



m=0 



where the numerical coefficients Vi m are related to the coefficients in Eq. (4.11) by expansion (4.18). The 
first five functions are 



Vi(t) 
Vs(t) 
V 4 (t) 
V 5 (t) 



45 2887 o 315 4 1463 fi 
I t + 4 j 1 

128 1920 128 1152 ' 

51 129 2 761 4 441 6 707 8 

+ t 2 1 4 + t 6 1 8 , 

128 32 64 32 128 ' 

3879 436931 2 264745 4 173209 fi 101395 „ 495271 in 



5120 35840 4608 1536 1024 ' ' 15360 

Integrating over z by parts and introducing polar coordinates, for the asymptotic part one finds 



,d-s-2 



+ r 2 ) l/2 Jo 



it/2 



d6sm d ~ 2 ecos 2m - s - 1 e. 



(4.24a) 
(4.24b) 
(4.24c) 
(4.24d) 

(4.25) 
(4.26) 

: / \ : 1 

1=1 m=0 

Due to the first beta function on the right this expression has a simple pole at s = — 1 with the residue 
-(«) -(-^fSi/vV g(m+l/2,(d-l)/2) 

p=0 m=0 w ^'^ ' 

where pd is defined in Eq. (2.37). Hence, the Laurent expansion of the asymptotic part near s = — 1 has 
the form 



M l roo 



The evaluation of the integrals leads to the following formula 

7 H, x s . vrs ^ ,wv^ T/ fd-s-1 l + s-d + 1 



s d-1 
B \ m , 

2' 2 



7 (as) 



f + 41o + 0(s + l), 



(4.28) 



with the finite part 

7 (as) _ 
^NKfi — 



2vrd 



p=0 



V'(p+l)-'0 o +V' ^+o -'0 " 7 '+o _2 



1 



/d-1 M \ I , 

y=l,d-Z=odd l=d+lj m=0 v 



2 

d Z-d 
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Table 2: Pole and finite parts of the Casimir energy for a single Neumann plate. 



d 


Si ^lNp 


,d-l 

SI IN f 


SI -^lATp 


Si IN f 




2 
3 
4 


t> 


-0.00874 
0.00213 
-0.000593 


5 


-0.000633 
0.000792 
-0.000398 


-0.00937 
0.00292 
-0.000194 


512tt(s+1) 
1 


512tt(s+1) 


180tt 2 (s+1) 
269 


180vr 2 (s+l) 
269 


262144vr 2 (s+l) 


262144tt 2 (s+1) 



where the second sum in the braces of the third line is present only for M > d + 1. Now adding the part 
coming from (4.22), for the Laurent expansion of the function Z^k{s) near s = — 1 one finds 



'NK 



{s) = 



7 {as) 
NK,— 

5 + 1 



y(as) y 



(1) 



(-l)+0(5 + l), 



(4.30) 



where the separate terms are defined by formulae (4.22), (4.27), (4.29). Using this result, for the vacuum 
energy in the case of Neumann scalar induced by a single plate at £ = £i in the region £ > £i one receives 



E 



w 4- R 



j \n (£u — -^lAfp 
where for the pole and finite contributions we have 



E 



(R) 
INp 



-l 



^(s + l)' 



E 



(R) 
IN} 



B d 
d-l 



4T 



(R) 
lAT/> 



7 (os) 7 i I i _ , 
NK,0 NK\ 1 ) 



(1) 



(4.31) 



(4.32) 



The numerical results corresponding to the vacuum energy (4.31) with separate pole and finite parts are 
presented in Table 2 for spatial dimensions d = 2, 3, 4. Comparing to the data from Table 1, we see that 
in all these cases the Neumann quantities for the RR region dominate the Dirichlet ones. 

By making use the formula given in Appendix B, we can obtain the vacuum energy in the RL region 
measured by an uniformly accelerated observer. In d = 3 spatial dimension for this energy one has 



W {R) ,£ N _ 9 
-kliVglSl) - 72 
Si 



0.00213 + 



1 



180vr 2 



1 



8 + 1 



57 J 



(4.33) 



It logarithmically depends on the normalization scale. 



4.2 Neumann vacuum energy in the RL region 

Similar to the case of the RR region, the Casimir energy for the Neumann scalar in the region between a 
single plate at £ = £i and the Rindler horizon corresponding to £ = is determined by the formula (see 
Sec. 5 below) 

b[ L n ) ^i) = ^iZ ni (s)\s=-i, (4.34) 
si 



with the function 



and the partial zeta function 



Z NI (s)= / dxx d - 2 c[ L J(s,x), (4.35) 
J o 

C$ (s, x) = i sin ^ J™ dz z~ s ^- In l' z (x). (4.36) 

For the analytic continuation of the function Z^j(s) to s = — 1 we note that for the derivative of the 
Bessel modified function one has the uniform asymptotic expansion 

I' z (x) = -^ ('' + za ) 1 V /^(4 /f )(x) ~ V Vl ® m , (4.37) 
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where the functions r/(x), t are defined in Eq. (2.16). Separating the contributions coming from the 
different factors in Eq. (4.37), we present the partial zeta function in the form 

(i L J(s,x) = ([ L N °\s,x) + (l L N 1 \s,x), (4.38) 

with 



C\s,x) = isin^ / dzz-± In- ^ +*) ^(s/,) (439) 
vr 2 J p dz V2vr x 

CS } ( S ,x) = isin^ r dzz- s ^\nl[ N \x). (4.40) 



For 1 < Res < 2 the z-integral in the expression for the function £^\s,x) is convergent in the limit 
p — > and the straightforward evaluation of the integral in this limit gives 



7r(l - s) 



Cj^s) = ~Z7i~~xB ( 1 - ^ ) sin^f + V- Ml) 



As in the case of RR region, the contribution of this term into the renormalized vacuum energy vanishes. 
Hence, we need to consider the quantity 

Zjvj(s) = Isiny y dxx d - 2 J dz z- s -^- In 4 N \x). (4.42) 

With the help of asymptotic expansion (4.37), it can be presented in the form 

Z NI ( S ) = Z^(s) + Z^ I (s), (4.43) 

where 

*£>(.) - L^£*^f^±£lM, ,4,4, 

zg, w _ ( ta ^ w _g_^). (445) 

Note that for M > d the expression (4.45) for Z^j(s) is finite at s = — 1 and we can directly put p = 0. 
After integrating by parts and introducing polar coordinates we find 

Z$\ N (-1) = ± /"drr"- 1 r /2 c» sin^^ [in (^e-W sin^ cose (r sin*)) - £ 

^ Jo Jo K ' l=1 (1 + r z ) ' 

(4.46) 

As regards the asymptotic part, in the expression for Z^f (s) we can directly put p = and after 
introducing polar coordinates and integrating one finds 

r,(as) , \ -5 . 7T5 v "v v "v ^ f d — S — 1 I + S — d + l\ ^ f S d — l\ , , 

4 a ?(*) = S - Y EE ^ ( — 5— . — 2 — ) B ( m -2—)- ^ 

1=1 m=0 v / \ / 

At s = — 1 this function has a simple pole with the residue 

4t-i = ("1)^,-1. ( 4 - 48 ) 
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with _ 1 defined by Eq. (4.27). Now, taking into account the contribution coming from (4.45), we 
have the following Laurent expansion near s = — 1 



(4.49) 



where 



7 {as) 
J NI,0 



1 



d-2p 



p=0 



m=0 



l/)(jp + 1) - tp 
d-l 



+ ij) ( m + 



d-l M \ l 

\J=l,d-Z=odd l=d+l m=0 



d l-d 



I d-l 
B[m + - 2 — 



(4.50) 



The vacuum energy for the Neumann scalar in the region < £ < £i is presented as a sum of the pole 
and finite terms 



p (L) p (L) 



with 



p (L) _ B d Z Nl]-l 
^lNp - t d-l, ' ' 



(i) 
IN} 



Ba 



y(as) y \ I! 

Z iVJ,0 + Z 7V/l _i J 



(4.51) 



(4.52) 



For spatial dimensions d = 2, 3, 4 the results of the corresponding numerical evaluations are reported in 
Table 2. For spatial dimension d = 3, the Casimir energy in the RL region measured by an uniformly 
accelerated observer is obtained with the help of formula (B.4): 



SI 



9 



-0.000792 - 



180vr 2 



8+1 



+ In 



(4.53) 



where g is the proper acceleration for the observer. 



4.3 Total Casimir energy for a single Neumann plate 

Summing the contributions from the RR and RL regions we obtain the total Casimir energy for a single 
Neumann plate: 



EMZi) = E{ N \Z 1 )+E[ L r }(£ 1 ). 
It can be divided into pole and finite parts, 

Ein(€i) = Ei Np + E 1N f, 

where 



tt\s + l) 



d-l 



Cl 



^NKfl + ^Nlfi + Z vM _i J + Z 7Vn _1 



/(as) 



(1) 



/(I), 



(4.54) 
(4.55) 

(4.56) 



In odd spatial dimensions the pole part vanishes due to the cancellation of corresponding RR and RL 
parts and the total Casimir energy is finite. In this case this energy can be presented in the form 



1 d-l , 
m + -, — - — ] + / drr 



2' 2 

Mi 



d-l 



w/2 



ln(-2rsin 2 W; cose (rsin^)K; cose (rsin^)-2^^ 
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Under the condition Mi > d/2 — 1 the integral in this formula is convergent. The numerical results for 
the finite part of the total vacuum energy are given in Table 2. In the physically most important case 
d = 3, the total Casimir energy for the Neumann scalar induced by a single plate is positive. For the 
total Casimir energy measured by an uniformly accelerated observer one obtains 

ElNt _ ~, (4.58) 

SI 

which is finite and independent on the normalization scale. As the energy (4.58) is positive, the corre- 
sponding vacuum forces tend to decelerate the plate. Here the situation is opposite compared to the case 
of the Dirichlet scalar. 



5 Casimir energy for two Neumann plates 

Consider the scalar vacuum in the region between two plates located at £ = £i, £ = £2 and with the 
Neumann boundary conditions on them: 

For the function </>(£) in expression (2.4) one has 

</>(£) = N iuJ (kt fc&) = C(H2)K^(kO - I^ikOKUkb). (5.2) 

From the boundary condition on the plate £ = £1 we obtain that the corresponding eigenfrequencies are 
solutions to the equation 

NUHi, fc£ 2 ) = iUk&Kiikti) - lUkti)K>Jk&) = 0. (5.3) 

We denote the positive roots to this equation by uj = k->Ar n (fc£i, /c£ 2 ), and will assume that they are 
arranged in the ascending order cu^n < WNn+i- F° r the vacuum energy in the region £1 < £ < £2 per unit 
surface of the plates on has 

1 f d d ^k 00 

EN= 2j (2vr) rf - 1 ^^ iVw(fcgl,fcg2) - (5 ' 4) 
We introduce the partial zeta function as 

00 

0v(s,fc£i,&£ 2 ) = J>^* ( 5 - 5 ) 

71=1 

In order to obtain the Casimir energy, one has to find the analytic continuation of (5.5) to s = — 1. 
Employing the Cauchy's theorem, zeta function (5.5) may be written under the form of contour integral 
on the complex plane, 

Gv(*,fc£i,fc£ 2 ) = ^ J dzz- a ^-]nNiMi,^2), (5.6) 

with the same contour C as in Eq. (3.7). The integral on the right of Eq. (5.6) can be presented in the 
form 

(N(s,Hi,k&) = C^\ S ,k^) + ^- j Ca dZZ ~ S §- z lnI -™( k & 



H — — ^ [ dzz s ^-ln 

2lTl ^ J C a oz 

a=+,— 



1 i'- m Mi)K[M2) 



+^-f c dzz- s ^-lnNUkZi,k&), (5.7) 
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where C + and C" are the upper and lower halves of the contour C except the parts coming from the 
small semicircle C p . The last integral over C p vanishes in the limit p — ► for Res < 2 and it can be 
omitted in the following consideration. After parametrizing the integrals over the imaginary axis we 
obtain the formula 

Gv(s,fc£i,fc6) = cS^O + ^sin^ J" ' dzz- s ^- z \nl'M2) 

1 . f°° ^ s d 

H — sin- / dzz s — 
vr 2 J p dz 

The last integral on the right of this formula is finite at s = — 1 and vanishes in the limits £i — ► and 
£2 — ► 00. It follows from here that the second term on the right corresponds to the zeta function ([^J (s, x) 
for the region on the left of a single plate located at £ = £2- The procedure for the analytic continuation 
of this function we have considered in the previous section. 

Taking into account (5.7), for the Neumann vacuum energy in the region £1 < £ < £2 one finds 

E N = + £7$ (&) + AENiZufr), (5.9) 

with the interference term 

r> f'OO f'QO 

A^at(6,6) = — / dkk d - 2 dzln 
n Jo Jo 

where we have integrated in parts the z-integral. From the inequality N' z (k^i, k^) < for £1 < £2 it 
follows that A£jv(£i,£2) < 0. In Fig. 4 we have presented the dependence of the interference part of the 
Casimir energy (5.10) on the ratio £1/^2 for d = 3. As seen from Figures 3 and 4, the interference parts 
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Figure 4: Interference part of the Casimir energy in the region between two Neumann plates, 
£2 _1 A£x>(£i, £2), as a function on the ratio £i/£2 for d = 3. 

of the Casimir energies for Dirichlet and Neumann scalars are numerically close to each other. This is 
a consequence of that the subintegrands in formulas (3.11) and (5.10) are numerically close. This can 
be also seen analytically by using the inequalities for the Bessel modified functions given in Ref. [17]. 
Interference term (5.10) is finite for all < £1 < £2, and diverges in the limit £1 — ► £2. In this limit the 
main contribution comes from the large values of z. Introducing a new integration variable x = k/z and 
using the uniform asymptotic expansions for the Bessel modified functions, for the interference part in 
the leading order one obtains the same result as in the Dirichlet case, (3.13). In particular, in the limit 
£i,£2 — ► 00 with fixed £2 — £1, for the vacuum energy measured by an uniformly accelerated observer with 
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the proper acceleration g — > and g£i — ► 1, one recovers the standard result for the Casimir plates on 
the Minkowski bulk. 

Adding to the energy in the region between the plates, given by Eq. (5.9), the energies coming from 
the regions £ < £i and £ > £2, we obtain the total Casimir energy, Ep 0t \ for two-plates geometry: 

E ( * ot) = £iaK£i) + £uv(6) + AEjv(6, 6), (5.11) 

where the interference part is given by formula (5.10). For d = 3 Neumann scalar the first two terms 
on the right of Eq. (5.11), corresponding to the single plates contributions, are positive and the third 
(interference) term is negative. For sufficiently close £1 and £2 the last term dominates and the total 
Casimir energy is negative. For the large separations between the plates the contributions from the 
single plate parts ar dominant and the total energy is positive. Hence, for some intermediate values of 
the separation the total Casimir energy for the Neumann scalar vanishes. 

As in the Dirichlet case, it can be seen that the energy (5.10) is related to the interaction forces 
between the plates by the relations 

= AU(a,6) = ^y o dkk ^ (5.12a) 

9 af (f r\ ^ f;n- Bd (°°dhh^ r, r ^(^i)(i + ^ 2 A 2 £ 2 2 ) 
^A^ (6 ,6) - -PW&,f 2 )--^y o dkk I dz r M 2)N>M.M 2 ) - (5 - 12b) 

Here p)yY int ), J = 1,2 are the vacuum effective pressures on the plate at £ = £j induced by the presence 
of the second plate. These quantities are investigated in Ref. [17]. 

6 Casimir energy for the electromagnetic field 

We now turn to the case of the electromagnetic field. We will assume that the plates are perfect conductors 
with the standard boundary conditions of vanishing of the normal component of the magnetic field and the 
tangential components of the electric field, evaluated at the local inertial frame in which the conductors 
are instantaneously at rest. As it has been shown in Ref. [15] for d = 3 and in Ref. [17] for an 
arbitrary spatial dimension (on the decomposition of the electromagnetic field in the Rindler coordinates 
see also Ref. [40]), the corresponding eigenfunctions for the vector potential may be resolved into 
one transverse magnetic (TM) mode and d — 2 transverse electric (TE) (with respect to the £ direction) 
modes A% a , cr = 0, 1, . . . , d — 2, a = (k, uj): 

At a = (-£^,-^,0,...,0^oa, <r = l, TMmode, (6.1) 
Ka = %<P*,«, <7 = 0,2,...,d-2, TE modes,. (6.2) 

Here the polarization vectors obey the following relations: 

e° a = el = 0, e afl e»,=-k 2 5 aa ,, e^ = 0. (6.3) 

From the perfect conductor boundary conditions one has the following conditions for the scalar fields 

(fiaa 1^=6 = V™ k=6 = °) cr = 0,l,...,d-2, (6.4) 

-or\t=*=-af\t=*=°- ( 6 - 5 ) 

As a result the TE/TM modes correspond to the Dirichlet/Neumann scalars and the Casimir energy for 
the electromagnetic field can be obtained from the scalar results by the formula 

E cm = (d- 2)E D + E N . (6.6) 
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The same relation takes place for the energies E\ e ^, E\J n , AE cm , Ecm ■ In particular, for the total 
electromagnetic Casimir energy in the geometry of a single plate from numerical results given in Tables 
1, 2 one has E em = 0.00160/^ in d = 3. Hence, in this case the electromagnetic forces tend to decelerate 
the plate. 

The total Casimir energy for two plates located at £ = £i and £ = £2 can be written in the form 

Et t} = £lcm(£l) + £lem(6) + AE em (^, (6.7) 

where AE em (^i,^2) < 0. Note that the quantity ^"^m' is a function on the ratio only. For a 

given £1 and for large values of this ratio the main contribution into E^^ comes from the first term on 
the right of Eq. (6.7). For £2/^1 — 1 1 the last term on the right of this formula dominates. From 
the numerical results given above it follows that in d = 3 for a given £1 the energy E^^ is positive for 
large values of ^2/^1 and is negative for ^2/^1 — 1 <C 1. Hence, for some intermediate value of £2/^1 this 
Casimir energy vanishes. 



7 Conclusion 

In quantum field theory the different unitary inequivalent representations of the commutation relations 
in general give rise to different pictures with different physical implications, in particular to different 
vacuum states. An interesting issue in the investigations of the Casimir effect is the dependence of the 
vacuum characteristics on the type of the vacuum. In this paper we have investigated the Casimir energies 
generated by a single and two parallel plates moving by uniform proper acceleration, assuming that the 
fields are prepared in the Fulling-Rindler vacuum state. The corresponding vacuum expectation values 
of the energy-momentum tensor were investigated in Refs. [15, 16] for the geometry of a single plate and 
in Ref. [17] in the case of two plates. Due to the well known surface divergencies in these expectation 
values, the total Casimir energy cannot be evaluated by direct integration of the vacuum energy density 
and needs an additional regularization. In this paper as a regularization method we employ the zeta 
function technique. We have considered the cases of scalar and electromagnetic fields in an arbitrary 
number of the spacetime dimensions. 

For the scalar case both Dirichlet and Neumann boundary conditions are investigated. In the case of a 
single plate geometry the right Rindler wedge is divided into two regions, referred as RR and RL regions. 
By using the Cauchiy's theorem on residues, we have constructed an integral representations for the zeta 
functions in both these regions, which are well suited for the analytic continuation. Subtracting and 
adding to the integrands leading terms of the corresponding uniform asymptotic expansions, we present 
the corresponding functions Z(s) as a sum of two parts. The first one is convergent at s = — 1 and can be 
evaluated numerically. In the second, asymptotic part the pole contributions are given explicitly in terms 
of beta function. As a consequence, the Casimir energies for separate RR and RL regions contain pole 
and finite contributions (see, for example, formulae (2.41) and (2.42) in the case of the Dirichlet scalar in 
the RR region). The remained pole term is a characteristic feature for the zeta function regularization 
method and has been found for many other cases of boundary geometries. The coefficient for this term 
is determined by the corresponding boundary coefficient in the heat kernel asymptotic expansion. For an 
infinitely thin plate taking RR and RL regions together, in odd spatial dimensions the pole parts cancel 
and the Casimir energy for the whole Rindler wedge is finite. Note that in this case the total Casimir 
energy can be directly evaluated by making use formulae (2.68) and (4.57) for the Dirichlet and Neumann 
scalar fields, respectively. The cancellation of the pole terms coming from oppositely oriented faces of 
infinitely thin smooth boundaries takes place in very many situations encountered in the literature. It is 
a simple consequence of the fact that the second fundamental forms are equal and opposite on the two 
faces of each boundary and, consequently, the value of the corresponding coefficient in the heat kernel 
expansion summed over the two faces of each boundary vanishes [21]. In even dimensions there is no 
such a cancellation. The numerical results for the separate pole and finite contributions to the Casimir 
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energy in RR and RL regions are summarized in Table 1 for the Dirichlet case and in Table 2 for the 
Neumann boundary condition. For the physically most important case d = 3 the total Casimir energy 
is negative for the Dirichlet scalar and positive for the Neumann scalar. This means that the vacuum 
forces tend to accelerate the Dirichlet plate and to decelerate the Neumann plate. 

In the case of two parallel plates configuration we have derived integral representations for both 
Dirichlet and Neumann zeta functions. The corresponding Casimir energies are presented as a sum of 
single plate parts and the interference term. The latter is determined by formula (3.11) for Dirichlet 
scalar and by formula (5.10) for Neumann scalar and is located in the region between the plates. It is 
always negative and is related to the corresponding interaction forces by the standard thermodynamical 
relations (see Eqs. (3.14), (5.12)). For large values of the separation between the plates, the total Casimir 
energy is dominated by the contributions coming from the single plate parts. For a given £i and small 
values of the separation the interference part is dominant and the total Casimir energy is negative. For 
d = 3 Dirichlet scalar this is the case for all values of the separation. However, for d = 3 Neumann scalar 
due to the positive single plate energies, the total Casimir energy for two plates is negative for small 
separations and positive for large separations. Consequently, this energy vanishes for some intermediate 
value of the ratio ^2/^1- Letting £2 - ► °o correspond to removing one of the plates and from the formulae 
for two plates we recover the Casimir energy for a single plate. The case d = 1 is considered separately 
in Appendix A. In this dimension the Casimir energies for Dirichlet and Neumann scalars are the same. 
They vanish for a single plate (point) geometry and are negative for the two plates case. Note that in 
d = 1 the problem is conformally related to the corresponding problem on background of the Minkowski 
spacetime. 

In Sec. 6 the case of the electromagnetic field is considered with the perfect conductor boundary con- 
ditions in the local inertial frame in which the boundaries are instantaneously at rest. The corresponding 
eigenmodes are superposition of TE and TM modes with Dirichlet and Neumann boundary conditions, 
respectively. The Casimir energies for the electromagnetic field can be derived from the corresponding 
scalar results making use formula (6.6). In particular, the total electromagnetic vacuum energy of a single 
plate in d = 3 is positive and the vacuum forces tend to decelerate the plate. For two plates geometry the 
situation for the electromagnetic field is similar to the Neumann scalar case. The total Casimir energy 
is negative for small values of the plates separation (the interference part dominates) and is positive for 
large separations (single plates parts dominate). 

In the main part of this paper we have considered the vacuum energy corresponding to the dimension- 
less coordinate r in Eq. (2.3). The corresponding eigenfrequencies are also dimensionless and there is no 
need to introduce an arbitrary mass scale in the definitions of the related zeta functions. However, this 
mass scale is necessary if we consider the vacuum energy measured by an uniformly accelerated observer. 
The relation between these energies is discussed in Appendix B, where we have shown that they are 
connected by formula (B.4). In this formula the logarithmic term with an arbitrary mass scale /x has to 
be viewed as a remainder of the renormalization process (for a discussion see [7, 21]). For infinitely thin 
boundaries in odd dimensions, in calculations of the total vacuum energy, including the parts from two 
sides of the boundary, these terms cancel and a unique result emerges. 
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A d = 1 case 



Let us consider the scalar vacuum in the region between two boundaries located at £ = £1 and £ = £2- 
The normalized eigenfunctions satisfying Dirichlet boundary conditions and the corresponding eigenfre- 
quencies are in the form [17] 

GXD ^^ti t 7TTI 

fn= ^^sin(w n ln(£ 2 /£)), u n = , n = 1,2, . . . . (A.l) 

For the corresponding vacuum energy one finds 

ED "2 n V" " ^(676) " "241n(6/6)' (A ' 2) 

where Cr( x ) is the Riemann zeta function. The renormalized Casimir energies E^J and E^j for a 
single boundary are obtained from (A.2) in the limits £1 — ► and £2 — ► 00, respectively, and both these 
quantities vanish. 

For the case of the Neumann boundary conditions the eigenfunctions have the form 

Pn = ^r= cos ^ M&/0) > n = 0, 1, 2, . . . (A.3) 
vvrn 



with the same eigenfrequencies as in Eq. (A.l). The corresponding Casimir energy is the same as in the 
Dirichlet case: En = Ejj- 



B Relation to the energy measured by an uniformly accelerated ob- 
server 

In this appendix we consider the relation of the energies evaluated above to the energy measured by an 
uniformly accelerated observer. The frequency u in Eq. (2.4) corresponds to the dimensionless coordinate 
r and, hence, is dimensionless. Proper time r 9 and the frequency lv 9 measured by an uniformly accelerated 
observer with the proper acceleration g and word line (x 1 ) 2 — t 2 = g~ 2 , are related to r and u by formulae 
Tg = t / g, ujg = Log (the features of the measurements for time, frequency, and length relative to a Rindler 
frame as compared to a Minkowski frame are discussed in Ref. [41]). The corresponding vacuum energy 
has the form 

1 r d d-i k - 

^ = 27(^J>- 

where oj gn = oj gn {x) are the eigenfrequencies measured by the uniformly accelerated observer. As above 
we introduce the partial zeta function related to these eigenfrequencies and the function Z g (s) by relations 

OO 

( g (s,x) = J2(" 9 n/vT S , Z g (s)= dxx d - 2 ( g (s,x). (B.2) 

n=l J ° 

Note that we have introduced an arbitrary scale fi with mass dimension, in order to keep the zeta function 
dimensionless for all s. Now the Casimir energy can be written as 

E g = ^±Z g (s)\ s= ^. (B.3) 

By taking into account the relation between the eigenfrequencies corresponding to r and r g coordinates, 
one obtains the relation between Z g (s) and the function Z(s) considered above: Z g (s) = (fi/g) s Z(s). 
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Laurent-expanding over s + 1 and using the corresponding expansion for the function Z(s) given above, 
we obtain the following relation between the Casimir energies: 

/ BhZ^ \ 

E g =glE + Mj,/g)-^\ , (B.4) 

where the energy E corresponds to the dimensionless coordinate r. In dependence of the boundary 
conditions and the region under consideration, here as a pair (E, Z^^) we have to take (E^) , Zp^ 

(E^p\ Z^pp_-^, (Eip, Zp 2 ^ _ 1 + Z^^), where -F = D,N,em. Note that in odd spatial dimensions for 
the total Casimir energy E\p g of a single plate the scale dependent parts from the RR and RL regions 
cancel and we obtain a scale independent energy. 
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